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Abstract 

The type-I intermittency route to (or out of) chaos is investigated within the Horizontal Visibility 
graph theory. For that purpose, we address the trajectories generated by unimodal maps close to 
^ \ an inverse tangent bifurcation and construct, according to the Horizontal Visibility algorithm, 

their associated graphs. We show how the alternation of laminar episodes and chaotic bursts has 
a fingerprint in the resulting graph structure. Accordingly, we derive a phenomenological theory 
. that predicts quantitative values of several network parameters. In particular, we predict that 

in 

OO ' the characteristic power law scaling of the mean length of laminar trend sizes is fully inherited 

in the variance of the graph degree distribution, in good agreement with the numerics. We also 

O ' 

. report numerical evidence on how the characteristic power-law scaling of the Lyapunov exponent 

as a function of the distance to the tangent bifurcation is inherited in the graph by an analogous 

^ ' scaling of the block entropy over the degree distribution. Furthermore, we are able to recast the 

H ■ 

' full set of HV graphs generated by intermittent dynamics into a renormalization group framework, 

where the fixed points of its graph-theoretical RG flow account for the different types of dynamics. 
We also establish that the nontrivial fixed point of this flow coincides with the tangency condition 
and that the corresponding invariant graph exhibit extremal entropic properties. 

PACS numbers: 05.45.Ac,05.45.Tp,89.75.Hc 
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I. INTRODUCTION 



One of the common transitions between regular and chaotic behavior is intermittency, 
the seemingly random alternation of long quasi-regular or laminar phases, so-called 
intermissions, and relatively short irregular or chaotic bursts. Intermittency is omnipresent 
in nonlinear science and has been weighed against comparable phenomena in nature, such 
as Belousov-Zhabotinski chemical reactions, Rayleigh-Benard instabilities, turbulence. 



etc. [l-4|. 



The study and characterization of the onset mechanisms and main statistical 



properties of intermittency was carried out already a long time ago; Pomeau and Man- 
neville jsl introduced a classification as types I-III for different kinds of intermittency. 
Subsequently, other types have been described and typified, such as on-off intermittency p , 
ring intermittency [3], etc. Our objective here is to generate networks from the time series 
associated with intermittency and look how this phenomenon translates into such a different 
setting, and then examine the manifestation of its properties in the new context. For 
definiteness we chose the case of type I intermittency as it occurs just preceding an (inverse) 
tangent bifurcation in nonlinear iterated maps, although the very same methodology can be 
extended to other situations. Specifically, we show how this phenomenon can be visualized 
through the graphs generated when the Horizontal Visibility (HV) algorithm IS 
applied to the trajectories of the universality class of unimodal maps, as represented by the 
quadratic logistic map. 



The idea o 



approaches 



mapping time series into graphs is actively developed at present via different 



g|-|l6l|. Amongst them, the HV approach offers a promising new method for 
performing time series analysis, most of all because it has been corroborated that the 
fundamental nature of rather different complex dynamical processes is inherited by the 
associated visibility graphs. As part of the effort of developing a mathematically-sound 
visibility graph theory of dynamical systems, in recent years the performance of the visibility 
method has been tested and found to be consistently capable in different circumstances, 
including the description of chaotic, fractal-stochastic, or dissipative processes, to cite 



some 



17 



18| . In every case the network counterpart of each particular kind of dynamics 
has been determined with precision, positioning the visibility algorithm as a well-defined 
method to analyze the dynamics of complex systems using graph-theoretical tools. In 
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the context of low- dimensional chaos, two main routes to chaos have been studied in the 
light of this technique. Specifically, the period-doubling bifurcation cascade (Feigenbaum 
scenario) and the quasiperiodic route have been analyzed through the HV formalism and 
two complete sets of graphs, called Feigenbaum and quasiperiodic graphs respectively, that 
encode the dynamics of their corresponding classes of iterated maps, have been introduced 
and characterized recently 19l422(|. The third well-known route to chaos present in 
low- dimensional dissipative systems is type-I intermittency (Pomeau-Manneville scenario), 
and in the present work we present the structural, scaling and entropic properties of the 
graphs obtained when the HV formalism is applied to this situation. 



In the following we first recall in section II the key aspects of type I intermittency. 
In section III we present the Horizontal Visibility algorithm and apply it to the study of 
trajectories generated by unimodal maps close to an inverse tangent bifurcation, where 
type-I intermittency takes place. A phenomenological derivation of the degree distribution 
P{k] e) of this kind of graphs is performed. We show that this single graph metric encodes 
the key scaling properties of type I intermittency, namely (i) the mean length {£) of the 
laminar episodes with e manifests in network space as a comparable scaling with the same 
variable of the second moment {k"^) of the degree distribution P{k,e), and (ii) the scaling 
of Lyapunov exponent A(e) is recovered in network space from the Shannon block entropies 
over P{k; e). 

Next, in section IV we recast the family of HV graphs generated by intermittent series into 
a graph-theoretical Renormalization Group (RG) framework and determine the RG fiows 
close to and at tangency. We show that there are two trivial fixed points akin to the high 
and low temperature fixed points in thermal phase transitions together with a nontrivial 
fixed point associated with the tangency condition. Finally, we also determine the extremal 
entropic properties of the RG fixed points as well as the entropy evolution along the RG 
fiows. 



II. TYPE-I INTERMITTENCY 



Type-I intermittency can be observed infinitely-many times in the logistic map 

xt+i = F{xt) = ixxt{l - xt), 0<a;<l, 0</i<4, 
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FIG. 1: Graphical illustration of how the Horizontal Visibility (HV) graph inherits in its structure 
the dynamics of the associated intermittent series. In the top of the figure we show a sample 
intermittent series generated by the logistic map close to /ic (e > 0), producing laminar regions 
(black) mixed with chaotic bursts (white). In the bottom we plot the associated HV graph. Laminar 
regions are mapped into nodes with a periodic backbone, whereas the actual pseudoperiodicity of 
the series is inherited in the graph by the existence of so called peak or interfacial nodes. Chaotic 
bursts are mapped into chaotic nodes, with a characteristic degree distribution (see the text). 



close to the control parameter values /i = /i^ at which windows of periodicity open with 
period T for values /i > /ioo — 3.569945672...., where /ioo is the accumulation point of the 
main period-doubling cascade that locates the first appearance of chaos when increasing 
/i from small values. It can be observed that at /is = 1 + y/S this map exhibits a cycle of 
period T = 3 with subsequent bifurcations. This is the most visible window of periodicity 
in the chaotic regime and the one in whose vicinity our simulations have been performed. 
The regular periodic orbits hold slightly above /ly but below fir the dynamics consists of 
laminar episodes interrupted by chaos, a phenomenon known as intermittency. In what 
follows we relabel fiT = /^c and define e = fic — fJ'- In the upper part of figure [1] we show a 
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sample type-I intermittent time series generated by the logistic map close to /is, showing 
alternation between laminar trends, represented by black dots in the series, and chaotic 
bursts, represented in turn by white dots. We note that the laminar phase is not actually 
periodic, but approaches a periodic behavior of period 3 and it is precisely this behavior 
which is close to periodicity what makes it easily distinguishable from the chaotic bursts. 



A. Basic properties of intermittent series generated by unimodal maps 

Under rather general circumstances, trajectories generated by canonical models evi- 
dencing type-I intermittency show power-law scaling in the mean length of laminar phases 
(£) ~ e~"' , where e, called the channel width of the Poincare section, is the distance between 



the local Poincare map and the diagonal 23|. The specific value of exponent 7 is typically 



associated to the reinjection mechanism and several exponents have been reported, although 



7 = 0.5 holds in a rather large set of situations 2^. In figure [2] we plot in log-log scales 
the size distribution of laminar phases P{i;e), derived numerically from time series of 10^ 
data generated through the logistic map, for different values of e, showing the characteristic 
asymmetric U-shape. In the bottom inset panel of the same figure, we plot the dependence 
of the mean length of laminar trends with e, showing the well studied scaling {£) ~ e~^/^. 
We finally note that the maximum length of the laminar trends ip, the extremum of the 
distribution P{i; e), scales also as £p a . 



On the other hand, the length of a chaotic burst is known to be unpredictable. In figure 
|3] we plot, in semi-log scales, the size distribution of chaotic bursts e) for the same 
time series as used in figure [2J This distribution has in turn an exponential decay which 
becomes fairly independent of e for sufficiently large lengths, a result that can be justified 
invoking the survival time of a random memoryless process. In the lower inset of figure 
|3] we numerically check that, as P{£b;e) ~ P{h), the mean length of the chaotic bursts 
remain constant independently of e, with an approximated value of {ib) ~ 15. 
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FIG. 2: Log-log plot of the length distribution of laminar phases P{i; e) in our system, for different 
values of e: (squares) e = 3 • 10"^, (up-triangles) e = 10^^, (down-triangles) e = 3 • 10"^, (right- 
triangles) e = 10^^, (left-triangles) e = 3 • 10~^, (diamonds) e = 10"'^, (circles) e = 3 • 10~^. 
The starting power law decay is saturated at certain lengths, obtaining the classical asymmetrical 
U-shaped curves. (Inset bottom panel) Log-log plot of the laminar phases mean length {£) as a 
function of e, yielding the expected scaling {£) ~ e"'^'^. (Inset upper panel) Log-log plot of the 
local maximum of the distribution Ip as a function of e, also yielding a scaling Ip ~ e"*^'^. 

III. TRANSFORMATION OF INTERMITTENT TIME SERIES INTO HORIZON- 
TAL VISIBILITY GRAPHS 

The Horizontal Visibility (HV) algorithm js, 9| assigns each datum Xj of a time series 
{xi}i=i^2,... to a node i in its associated HV graph (HVg), where i and j are two connected 
nodes if Xi,Xj > x„ for all n such that i < n < j. Structural properties of a time series 
are inherited by its HVg (see the appendix for some specific properties relevant to the 
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FIG. 3: Semi-log plot of the size distribution of chaotic bursts in the system, showing a clear 
exponential decay which is fairly independent of e. This behavior, which can be understood as a 
result of the survival time of an random uncorrelated process, suggests that the relaminarization 
process occurs at unpredictable, uncorrelated times. (Inset bottom panel) Mean length of the 
chaotic bursts (/;,) for different values of e, showing a constant value (/;,) 15. (Inset upper 
panel) Semi-log plot of the degree distribution associated to the nodes belonging to chaotic bursts 
Pb{k), for different values of e. The exponential decay, characteristic of an HV graph associated 
to a chaotic series, has been already found in several other maps P, The results seem to 
be independent of e, in good agreement with the fact that the relaminarization process occurs at 
unpredictable times, what suggests Pb{k) = Aex.p{—ck). 

intermittent structure). In the bottom part of figured] we show the HV graph of the asso- 
ciated intermittent series, which consists of several repetitions of a 3-node motif (periodic 
backbone) hnked to the first node of the subsequent laminar trend, interwoven with groups 
of nodes irregularly (chaotically) connected amongst them. We observe that the motif repe- 
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titions in the graph correspond to the laminar regions in the trajectory (pseudoperiodic data 
with pseudoperiod 3) and the chaotically connected groups correspond to the chaotic bursts 
in the trajectory. As laminar trends are indeed pseudoperiodic in the sense that they can be 
decomposed as a periodic signal and a drift (see appendix), this pseudoperiodicity expresses 
in the graph structure by allowing a node for each period-3 motif to be connected to the 
first node in the next laminar region (the so called peak or interfacial node), as the values 
of the time series in the chaotic bursts are always smaller than those in the former laminar 
trend. The sequence of degrees is of the form 2 — 3 — 6 for laminar trends and loses this 
pattern in the chaotic burst. At odds with standard approaches, for which the distinction 
between laminar and chaotic phases is somewhat ambiguous, in this work we take advantage 
of this characteristic pattern as the criterion to numerically distinguish between both phases. 



The assembly of repeated 3-node motifs separated by sets of nodes with chaotic links 
inherited by the HVg from the laminar trends and chaotic bursts in the intermittent series, 
leaves also a characteristic footprint in its degree distribution P{k;e) (see figure Hj). This 
can be seen when P{k; e) is compared with that for fully chaotic motion (see for instance 



171]). In the former the connectivity of the nodes in the 3-node motifs is overrepresented 



than in the latter, since their relative frequencies are proportional to the length of the 
laminar episodes. Also, nodes with large degree (peak or interfacial nodes), increasingly 
large as e decreases, emerge due to reinjections after chaotic bursts because these have 
visibility over the nodes from laminar phases. The evidence collected leads us to express 
the total degree distribution P{k; e) as composed of three contributions that originate from 
three different types of nodes. Namely: 



i) Laminar: the contribution from the laminar phases, Pi{k;e), consists of a discrete 
set of degrees that correspond to the background periodic behavior. By construction (see 
appendix) a periodic series with a superimposed drift generates a graph where the nodes 
have, for the particular case in which T = 3, a degree belonging to the set {2, 3, 6}, therefore 
Pi{k] e) = Pi{k) = 1/3 when k = 2,3, 6, and zero otherwise. 

ii) Chaotic: the contribution from the chaotic bursts, P},{k;e), which, according to 
previous works j^, has an exponential decay (see the upper inset panel of figure [3] 
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FIG. 4: Log-log plot of the total degree distribution P{k; e) of 7 HV graphs mapped from trajec- 
tories of the logistic map in the vicinity of the window of period 3 (e > 0). Squares: e = 3 • 10~^, 
Up-triangles: e = 10^^, Down-triangles: e = 3 • 10~^, right-triangles: e = 10^^, left-triangles: 
e = 3 • 10^'', diamonds: e = 10"'', circles: e = 3 • 10"*^. Note that the tail of each distribution, 
associated to the peak nodes, scales with e'^'^, while this scaling is absent in the tail of the peak 
degree distribution (see the main panel of figure [5]) . This scaling is reminiscent of the contribution 
of the amount of peak nodes present in each series, whose abundance scales with e'^'^ (see the inset 
upper panel of this figure). (Inset upper panel) Log-log plot of the fraction of nodes associated with 
chaotic bursts, showing a scaling /fe(e) ~ 54 •e'^'^ (circles), and fractions of peak (interfacial) nodes, 
showing a scaling /p(e) ~ 3.2 • e'^-^ (squares). (Inset bottom panel) Log-log plot of the variance of 
the total degree distribution o"^ = (A;^) — {k)'^ as a function of e, obtained from the same graphs 
whose degree distribution is plotted in the main panel. A power law scaling of the form o"^ ~ g-o.5 
is found, which is the graph analogue of the well-known laminar phase scaling reported in the inset 
bottom panel of figure [2j 
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FIG. 5: Log-log plot of the degree distribution of the so called peak nodes, obtained from the 
HV graphs at different values of e. The asymmetrical U-shape is qualitatively similar to the size 
distribution of laminar phases (see the text). (Inset bottom panel) Log- log plot of the mean degree 
of peak nodes, as a function of e, calculated from the distributions plotted in the main panel of 
the figure. The scaling coincides with {I) /2> (see the text). (Inset upper panel) Log-log plot of the 
local maximum of the distribution kp as a function of e showing a scaling kp ~ e~'^'^. 

for numerical evidence). Moreover, as argued in the previous section, since the mean 
size of chaotic bursts is independent of e (as we can see in the inset bottom panel of 
figure [3]), the contribution Pb{k;e) is thus also independent of e, and in general reads 
Pb{k] e) = Pb{k) ~ exp(— c/c), where the specific value of c depends on the chaos dimension- 



ality [17 1 . 



iii) Peak: the contribution from the interface between the chaotic and laminar phases, 
Pp{k;e) (see figure [5]), arises from the peak nodes, with a very large degree that is 
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approximately proportional to the size of the laminar phase. Roughly, all of the peak 
nodes inherit a degree based on their visibility of the previous laminar phase, and therefore 
we a priori may assume that Pp{k;e) and P{i/3;e) (see figures [2] and E]). The factor of 
1/3 arises because, for the periodic window studied, k ~ i/3 (or k ~ i/T if the study 
focuses on intermittency close to a periodic window of period T), as we can observe in 
figure [H In particular, the maximum laminar size and peak node degree scale similarly, 
ip/3 ^ kp e~°-^ (see the upper panels of figures [2] and [5]). Note however that the 
connectivity k of the peak nodes is not straightforwardly distributed as i, since the actual 
value of the degree of the peak node assigned by the HV is not necessarily equal to 
the size of the full laminar phase, but varies in somewhat according to the actual value 
(position) of the series datum at reinjection (compare figures O and [2]). That is to say, 
visibility of a peak node may be as large as the preceding laminar trend, but it can 
also be smaller if reinjection takes place below tangency, or larger, allowing full visibility 
of the preceding trend and part of the previous one, if reinjection takes place above tangency. 

Nevertheless, in the following we will argue, and have numerically checked, that Pp{k; e) 
and P(£/3;e) have, up to first order, similar first and second moments, concretely: (a) 
The mean degree of peak nodes is approximately equal to 1/3-rd the mean value of the 
laminar phase size. This is due to the fact that the degree of peak nodes statistically 
self-averages over laminar phase sizes, as the variability in reinjection is symmetrical with 
respect to tangency. And (b) the second moment of Pp{k; e) is associated with both the 
second moment of P{i; e) and the variance of the reinjection distribution. However, the 
contribution of the latter is typically much smaller than the former. Indeed, whereas fluctu- 
ations in the reinjection position tend to decrease as e decreases, fluctuations in the size of 
laminar trends tend to increase as e decreases (this has been confirmed numerically) . There- 
fore, the leading contribution comes from the variability of laminar sizes for small values of e. 

The aforementioned phenomenology lead us to formally write down P{k;e) as 

P{k; e) = Me) ■ Pi{k) + /,(e) ■ P,{k) + /,(e) ■ Pp{k- e), (2) 

where fi{e),fb{e) and /p(e) are the fractions of nodes in the graph that correspond to 
laminar, chaotic and peak regions respectively. In order to derive these fractions, we rely 
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on three restrictions, namely: 



(i) Normalization: which trivially implies fi{e) + fb{e) + fp{e) = 1, 

(ii) Bounded degree: it has been proved jo], [lol, 3 that aperiodic series generate an HVg 
with constant mean degree (k) = 4, the upper bound value for HV graphs. This restriction 
implies 

Me) ■ J2 kPiik) + Me) ■ J] kP,{k) + /,(e) • J] kP,{k; e) = 4. (3) 

k k k 

Note that in the latter expression, the first sum is trivially Ylk^Pii^) = 11/3, the 
second sum corresponds to the graph associated with a chaotic series, and is directly 
^i^kPf,{k) = 4 [lol, Q, while the third sum yields, due to the aforementioned arguments, 
E,fcP,(A;;e)^(£)/3. 

(iii) After each chaotic burst a peak node emerges and anticipates the next laminar region, 
what implies fb = fp- (^b), where (£b) is the mean size of a chaotic burst and has been 
argued to be e- independent. 

After a little algebra, (i), (ii) and (iii) along with the rest of the arguments yield the 
prediction (4) = 11, A = H ■ M // = 1 - 12 ■ /p and fp = {£)-\ 

In order to compare our phenomenological prediction with the numerics, we recall that 



under rather general conditions {£) ~ 0.2 ■ e~°'^ 26|], which indeed coincides with the 
scaling plotted in the lower inset panel of figure [2l The predicted values for the fractions 
are therefore fp = 5 ■ e^'^ , fb = 55 ■ e"'^, and fi = 1 — 60 ■ e°'^. These can be compared 
with the results of numerical simulations, shown in the upper inset panel of figure El 
for which the best fit are fp ^ 3.2 ■ e°-^, ^ 54 ■ e°-^, fi^l — 57.2 ■ e°-^, and addition- 
ally (4) ~ 15 (bottom inset panel of figure [3]), on fairly good agreement with our prediction. 



A. Variance = (fc^) — (A;)^: graph analogue of {£) 



While there can be no equivalence between {£) anc 
latter is fixed to be {k) = 4 for an aperiodic regime 18| 



(k) in the intermittent graphs as the 
, a relationship may hold between {£) 
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and higher moments of P{k; e). Note that there is an increasing dispersion of the values of 
k from its mean (k) in the degree distributions of the graphs as e — )■ O"'', that we show below 
is related to the distribution connectivity Pp{k; e). If we measure this dispersion by means 
of the variance of the total degree distribution cr| = (/c^) — (A;)^ we recover numerically the 
g-0.5 scaling^ as shown in the bottom inset panel of figure HI 



This scaling is also a prediction of our phenomenological theory: 

4 = {k') - {kf 

= fi{e) J2 k'Piik) + Me) J] k'P.ik) + /,(e) J] PP,{k; e) - {k)\ (4) 

k k k 

Proceeding as before, the first sum is '^^k'^Pi^k) = 49/3 and the second sum '^^k'^Pb^k) 
is finite and e-independent. To determine the third sum we recall that the variance the 
distributions Pp{k;e) and P{i]€) are equal up to first order in (the variance associ- 
ated with the reinjection probability is always a smaller quantity), yielding k'^Ppik] e) ~ 
^^£^P(£;e) ~ e~^, where the higher-order terms take into account the fluctuations associ- 
ated to the reinjection variability. Collecting these results we straightforwardly find 

^2 ^-0.5 /c\ 

^fc ~ e , (5) 

for small values of e. This is a main quantitative result linking the key property of intermit- 
tent time series with its counterpart in the corresponding HV graphs. 

To round off, we here provide as an ansatz a concrete expression for Pp{k\ e): 

where kp = ae~^'^. This ansatz complies with all the properties obtained in our phenomeno- 
logical theory (although note that we do not need to use it to derive several observables as 
the fraction of nodes or the variance al), and provide a closed expression for the total degree 
distribution. This closed expression allows us to calculate the concrete degree distribution 
at tangency. This consists of an infinite-size laminar phase and a 'phantom' of kp with di- 
verging degree and vanishing probability of occurrence. Incidentally, note that in this limit 
case we also recover {k) = ^ kP{k; e — )■ 0) = 4, as expected. 
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FIG. 6: Log-log plot of the block entropies hn constructed from degree distributions of n-sequence 
of connectivities in the HV graphs as a function of e: n = 1 (squares), n = 2 (up-triangles) , n = 3 
(down-triangles), n = 4 (right-triangles). A scaling of the form hn ~ e"^"'^ reminiscent of the well 
known scaling of the Lyapunov exponent A ~ is found, albeit with a different, n dependent 
exponent a{n). (Inset panel) Log-log plot of the convergence of a{n) to the exponent associated to 
the Lyapunov exponent, as a function of n. A relation of the form [0.5 — a(n)] ~ yi~0.19 found, 
whose extrapolation suggests lim„_>oo a{n) = 0.5 



B. Scaling of Lyapunov exponent: Block entropies /i„ 



In the previous section we have studied how the scaling of (/) was inherited in P{k) by 
kp and al, but there exists another well known scaling relation in the type-I intermittency 
route involving the Lyapunov exponent 

71-1 



lim -V ln|F'(a;i)| (7) 

rt— >oo n — ^ 



A 

n 



1=0 

^0.5 



of the trajectories [5|,|25[, which reads A ~ e as e — 0. We recall here that the Pesin identity 
relates the positive Lyapunov exponents of a chaotic dynamics with the Kolmogorov- Sinai 



entropy of the system. In recent works 



19 



22| . such relation has been investigated 



in the visibility graph framework, through the definition of a graph theoretical entropy, a 
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Shannon-like entropy over the degree distribution 

^i = -$^P(A;)logP(A;). (8) 

k 

In figure [6] we plot in log-log the values of hi (solid squares) as a function of the channel 
width e. A power law scaling is recovered, albeit with a different scaling exponent a < 
0.5. Notice however that hi is only a proxy of the Kolmogorov-Sinai entropy and thus a 
comparison with the Lyapunov exponent is only approximate. Interestingly, hi is indeed 
the graph theoretical version of a size-1 block entropy over the degree distribution. Since 
the Kolmogorov-Sinai entropy of a map can be recovered as the asymptotic limit of block 



entropies 



28| s(n) = — ;i J2xi a;„P(^i5 ■■■y^n) logp(xi, ...,x„), wc take advantage of this fact 



to define a set of graph-theoretical block entropies 



hr 



-- V P{ki,...,K)iogP{ki,...,K). (9) 



n 



The e-dependence of these block entropies are also plotted in figure for different block sizes 
n. We observe in every case a power law scaling /i„ ~ . In the inset of the same figure 
we show how the exponent a{n) converges with n to the exponent found for the scaling of 
the Lyapunov exponent, suggesting 

lim hr, = A. (10) 

n—>oo 

We remark at this point that, whereas the entropy is a magnitude defined in the graph, 
the Lyapunov exponent is only defined in the system. Therefore, in rigor the Pesin identity 
cannot be used here as the explanation for equation [TOl as we are mixing properties defined 
in two different contexts. However, the strong numerical evidence in favor of a Pesin-like 
identity between the map's Lyapunov exponent and the graph's block entropy sug gest that 



a graph analogue of the Lyapunov exponent may be defined in network space 



221. 



IV. GRAPH-THEORETICAL RENORMALIZATION GROUP ANALYSIS 

Once the characteristic scalings of type I intermittency have been recovered in the 
topology of the asssociated graphs, let us provide a wider picture of the phenomenon. The 
overall properties of intermittency graphs can be framed in the context of a Renormalization 



Group (RG) transformation by following the procedure of previous studies |l9l-l2l|. We 
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FIG. 7: Illustration of the renormalization operator TZ applied on the HV graph at e = 0. This 
graph renormalizes, after two iterations of TZ, into an HV graph Gc which is itself (i) invariant 
under TZ, and (ii) unstable under perturbations in e, thus constituting a nontrivial (saddle) fixed 
point of the graph-theoretical RG flow. 

define an RG transformation TZ on an HV graph G as the coarse-graining of every couple of 
adjacent nodes where at least one of them has degree k = 2, into a block node that inherits 
the links of the previous two nodes. Iterating TZ we can trace the RG flows of intermittent 
graphs G{e). Results include: 



(i) When e < (yU > fic) trajectories are periodic and every HV graph trivially renormalizes 



towards the so called chain graph Gq, an infinite chain with k = 2 for all nodes |19l. I20[|. Go 
is invariant under renormalization TZ{Go} = Go, and indeed constitutes a trivial (attractive) 
fixed point of the RG flow, 7^("HG'(e < 0)} = Gq. 



(ii) When e > (/i < fic) repeated RG transformations eliminate progressively the links 
in the graph associated with correlations in the time series, leading ultimately to the 
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6>o e<o 

FIG. 8: Graph-theoretical RG flow of intermittent graphs. Graphs renormalize to Gq for e < 
and to Grand for e > 0, the trivial (attractive) fixed points of the flow. The graph obtained at 
e = renormalizes into Gc, which invariant under the RG transformation and unstable under 
perturbations in e, constituting the nontrivial graph-theoretical fixed point of the dynamics, whose 
degree distribution coincides with Eq. ([T2|) for e = 0. 

HV graph that corresponds to a random time series. The links that stem from temporal 
correlated data connect primarily laminar nodes, whereas the links between either burst 
and peak nodes originate from uncorrelated segments of the time series. If the laminar 
episodes are eliminated from the time series, the burst and reinjection data values form 
a new time series that upon renormalization leads to the random time series. We have 



7^("){G(e > 0)} = Grand , where G^and is t 
uncorrelated process with known graph properties 
(attractive) fixed point of the RG flow. 



le HVg associated with a random 



20|. This constitutes the second 



(iii) At e = (/i = fic) the HV graph generated by trajectories at tangency con- 
verges after only two steps of the RG transformation to a nontrivial flxed point 
TZ^{G{e = 0)} = Gc = TZ{Gc} and remains invariant under TZ afterwards. This feature can 
be demonstrated by explicit application of 71 upon G(e = 0) (see flgure[7]for a graphical 
illustration of this process). The flxed-point graph Gc is the HVg of a monotonically de- 
creasing time series bounded at inflnity by a large value, Gc is unstable under perturbations 
in e and it is thus technically a saddle-point of the RG flow, attractive along the critical 
manifold (spanned by G(e = 0) and its replicas within other periodic windows of period T). 
The RG flow diagram is shown in Figure [HI 
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From visual inspection of figure [7] the degree distribution for G{e = 0) is by construction 



P{k;0) 



lim < 

Af-s>oo 



N/3 
N 



if = 2,3, 



N 



if A; = 6, 



(11) 



i if 



N/3, 



and zero otherwise, where is the number of nodes. Note that this expression coincides 
with the one predicted from our phenomeno logical theory when e — ?■ 0. Its first moment 
is (k) = 4 while its second moment (fc^) diverges, as expected from the previous result 



-0.5 



. By construction the degree distribution of Gc is 



N-2 
N 



if = 3, 



lim < 



i if A; = 2, 



(12) 



i if k 



N 



and zero otherwise. The mean degree is again {k) = 4 and its second moment also diverges. 



Let us return now, within our RG treatment, to the concept of graph-theoretical entropies 
hn- We first look at the dependence of hn on e. This is shown when e > in figure [6] in 
logarithmic scales and observe that these functions are power laws that reach in every case 
a minimum at tangency. For concreteness we focus on hi, for which hi{e = 0) = log 3. 
This minimum value of hi is retained for all e < (but with |e| below the period-doubling 
bifurcations that take place within the window of period three). Hence, entropy reaches a 
global minimum for the HV graph at tangency e = 0. Next we enquire about the effect 
of the RG transformations on hi. The entropy at the nontrivial fixed point vanishes, as 
hi [Pgc(^)] ~^ when N oo, that is, the RG reduces hi when e = 0. Also, the RG 



transformations increases hi when e >_0 
it when e < (since hi [Pcoik)] = 



(as hi[PG,,^^ik)] = log(27/4) (l9|,|20|) and it reduces 
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20|). When e > the renormalization process of 
removal at each stage of all nodes with k = 2 eliminates node-node correlations (temporal 
correlations of the intermittent dynamics) and leads to a limiting renormalized system that 
consists only of a collection of uncorrelated variables, generating an irreversible fiow along 
which the entropy grows. On the other hand, when e < renormalization increments the 
fraction of nodes with degree /c = 2 at each stage driving the graph structure towards the 
simple chain Gq and thus decreases its entropy to its minimum value. Thus we observe the 
familiar picture of the RG treatment of a model phase transition, two trivial fixed points 

18 



that represent disordered and ordered, or high and low temperature, phases, and a nontrivial 
fixed point with scale-invariant properties that represents the critical point. There is only 
one relevant variable, e = ^c~fJ'j that is necessary to vanish to enable the RG transformation 
to access the nontrivial fixed point. The property that is seldom observed [29] is that an 
entropy functional, in the present case hi[P{k; e)], varies monotonously along the RG fiows 
and is extremal at the fixed points. A salient feature of the HV studies of the routes to 



chaos in low- dimensional nonlinear iterated maps, period doubling 
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20], quasiperiodicity 



2l|, and intermittency as presented here, is the demonstration that the entropy functional 
hi[P{k)] behaves as mentioned and attains extremal (maxima, minima or saddle-point) 
values at the RG fixed points. 

V. SUMMARY 

We have demonstrated the capability of the HV algorithm for transforming into network 
language the properties of the route to chaos via intermittency of type I as it occurs 
in unimodal one-dimensional iterated maps. The outcome is a novel type of network 
architecture composed by nodes of three types with characteristic connectivities and that 
are the building blocks with which intermittency is expressed recursively via concatenation 
in network space. These node types arise from the laminar trends, the chaotic bursts, and 
the interfacial positions related to reinjection into the channel. Their relative numbers 
as a function of the channel width e and their contributions to the degree distribution 
P{k) were determined from a phenomenological theory, with accurate results in agreement 
with numerical simulations. We have shown that the characteristic ingredients of type-I 
intermittency are inherited by the graph as functionals of the latter degree distribution, 
namely (i) the graph analogue of the laminar phase mean length, which evidences the 
well known scaling {£) ~ e^*^'^, was identified to be the variance of the total degree 
distribution, and (ii) the Shannon block entropy over the degree distribution appears as 
the graph analogue of the Lyapunov exponent A, as the characteristic scaling A ~ g-o-5 jg 
asymptotically recovered in network space through graph-theoretical block entropies hn- 
We note at this point that Pesin identity suggests that a graph theoretical analogue to 



the Lyapunov exponents could be defined in the network context 



221 ]. We also would like 



to highlight that while we have focused, for definiteness, in the logistic map close to the 
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periodic window of period 3, results can be trivially extended to the intermittent dynamics 
close to any given periodic window or any chaotic map undergoing an inverse tangent 
bifurcation. 



Significantly, the HV formalism leads to analytical expressions for the degree distri- 
bution near and at tangency. The scaling properties of the intermittent networks can be 
determined in terms of the same RG transformation employed with success on the HV 



graphs obtained for the period doubling and quasiperiodicity routes to chaos |19l-l2l[|. The 
graph-theoretical RG fixed points capture the features of the dynamics above, below and 
at the tangent bifurcation. Finally, the optimization of a graph entropy introduced via the 
degree distribution reproduces the RG flows and fixed points. 



In conclusion, the transition to chaos via type-I intermittency, as exemplified in uni- 
modal maps near and at an inverse tangent bifurcation, has been fully described within the 
Horizontal Visibility theory. This technique may be extended to the study of other types 
of intermittency and may be useful for the analysis and interpretation of time series with 
sporadic features of diverse origin. 
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VI. APPENDIX 



Some properties of HV graphs which are relevant to the characteristics of intermittent 
series include: 



i) Drift: Let {xi}i=i^2,... be a monotonically increasing/decreasing series and define 
its drift as {di}i-i^2,... = {^i+i — 2;i}i=i,2,...- Every node i of its HVg is connected only 
to the previous node i — 1 and to the following node i + 1, as it is always true that 
Xi+i ^ Xj ^ Xi-i or Xj+i ^ a;j ^ Therefore every node has connectivity k = 2 and the 
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degree distribution is P{k = 2) = 1. 



ii) Periodicity: Let {xj}j=i^2,... be a periodic series of period T so that Xj+y = x,. 
Let Xmax be max{xfc}fc=i,...,i+T- Its HVg consists of a repeated periodic motif of T nodes 



from node i^ax corresponding to Xmax to node imax + T corresponding to x^ax+T = x^ax [18 1. 



iii) Periodicity with drift: Let {xf}i=i^2,... be a drifted periodic series with drifted period 
T so that xf^rp = xf + di. If the series satisfies | xfj.^ — xf \>\ di — rfj+i | for all i, the HV 
graphs associated to {xf} and to {xj} in ii) are identical. 

iv) Chaoticity: Let {xj}i=i^2,... be a chaotic series. The degree distribution of its HVg 
has an exponential tail P{k) ~ exp(— Afc) with A ^ ln(3/2), the specific value of A depends 
on the chaotic process from which the series has been extracted from [l^] . 

Also note that the initial period-3 cycle behaves as X2 = /XcXi(l — xi), xs = fJ.cX2{^ — X2), 
xi = ^0X3(1 — X3), and in the chaotic vicinity of this orbit we have 

X2 = (/ic — — Xi) = — Xi) — e[xi(l — Xi)] = X2 — efi~^X2, (13) 

^3 = il^c - e)a;2(l - eAtc^^)[l - 2:2(1 - e/Uc ^)] ~ 2:3 - e/i^^(2x3 - fXcxj), (14) 
xf ~ Xi — e/i~"'^(xi — 4yUcX3 + 2/icX3 + 2/1^x3X3 — jJ^^xl) (15) 

that can be seen to be a perturbed periodic orbit. If we consider the series {xf }j=i^2,... to 
be periodic with a drift, xf_^_rp = xf + di, then according to condition (iii) above the periodic 
HVg remains invariant for pseudoperiodic orbits {x^}j=i 2,.... 
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